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Abstract 

An approximation to Einstein's field equations in Arnowitt-Deser-Misner (ADM) canonical for- 
malism is presented which corresponds to the magneto-hydrodynamics (MHD) approximation in 
electrodynamics. It results in coupled elliptic equations which represent the maximum of elliptic- 
type structure of Einstein's theory and naturally generalizes previous conformal-flat truncations of 
the theory. The Hamiltonian, in this approximation, is identical with the non-dissipative part of the 
Einsteinian one through the third post-Newtonian order. The proposed scheme, where stationary 
spacetimes are exactly reproduced, should be useful to construct realistic initial data for general 
relativistic simulations as well as to model astrophysical scenarios, where gravitational radiation 
reaction can be neglected. 
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I. INTRODUCTION 



In electrodynamics, to model quasi- stationary scenarios accurately, it is customary to 
drop the Maxwell displacement current in Maxwell- Ampere law of Maxwell's equations. This 
approximation is often referred to as the magneto-hydrodynamics (MHD) approximation as 
it is mainly implemented in MHD computations dealing, mostly, with astrophysical plasmas. 
A strictly analogous approximation to Einstein's general relativity, should prove very useful 
to model non-stationary situations, involving relativistic self gravitating systems. 

In Einstein's theory of gravity, the conformal-flat approximations played a major role 
within numerical relativity as they resulted in simple elliptic equations Q]. However, in 
the past, it was pointed out that these approximations are rather crude for highly non- 
spherical objects like rotating disks of dust or compact binary systems Therefore, 
in the following, we shall propose an elliptic-type approximation which goes far beyond 
conformal-flat schemes. Whereas the conformal-flat approximations are not exact under 
stationarity conditions, in the new approach, Einstein stationarity is reproduced exactly. 
In a post-Newtonian setting, the new approach, in general, does not reproduce all even- 
order post-Newtonian terms (of course, the odd terms are not part of the new approach) 
beyond the second post-Newtonian order. However, at the third post-Newtonian order 
of approximation, it still gives the equations of motion, derivable from the full Einstein 
theory 0, 0] . We observe that the very recently proposed gravito-anelastic approximation is 
quite different from our minimal no-radiation approximation, especially for the choice of the 
dropped and the kept time derivatives in the evolution equations Q]. Our approximation 
also results in elliptic equations, which should be solvable employing LORENE, a highly 
efficient numerical library, based on spectral methods, developed by the Meudon group [7]. 



II. EINSTEIN'S THEORY IN THE ADM FORMALISM 

In the ADM formulation of general relativity, the spacetime line element in the (3 + 1) 
decomposed form is given by 

ds 2 = -a 2 c 2 dt 2 + -f i:j (dx l + i3 i cdt){dx j + p j cdt), ( i, j = 1, 2, 3 ) , (1) 

where a is the lapse function, f3 l the shift vector, 7^ the induced metric on a three- 
dimensional spatial slice £(£), parametrized by the time coordinate t, and c is the velocity of 
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light. The 3-metric and its canonical conjugate yf-Q vr 4 - 7 , which is a contravariant symmetric 
tensor density of weight +1, satisfy the Hamiltonian and momentum constraints 

7 1/2 R = - - 2 ) + ^7 1/2 « 2 T 00 , (2) 



-2<. + vr% M = i^ 7 1 /^0 ) (3) 

where R is the curvature scalar of £(£), 7 the determinant of 7^, 7rj = 7jfc7r ifc , n = 7r|. 
T 00 and are the components of the unspecified 4-dimensional stress-energy tensor for 
the matter, T^ v . The canonical conjugate ^ is related to Kij, the extrinsic curvature of 
E(i), by 7T*- 3 = — 7 1 / 2 (7*'7 jm — 7 lJ 7' m )^ m , where 7*' is the inverse metric of 7^'. In the 
above equations, a partial derivative is denoted by , and G is the Newtonian gravitational 
constant. We note that the more familiar form for the left hand side of Eq. (j3*|), namely 
irlj, where | stands for the 3-dimensional covariant derivative, is expanded as 7rf 3 - — \^ kl 1ki,ii 
using the fact that is a mixed tensor density of weight +1. The 3-metric and its canonical 
conjugate evolve in accordance with the following evolution equations [s| 

7r 4 % = ~a-f 1/2 (2R ij - 7 ij R) + ia7" 1/2 T y (2C< ~ " a7~ 1/2 (27r im < - ttt^) 

+^ a 7 1/2 7 l '7 jm T im (4) 

and 

7y,o = a7 _1/2 (27Tj i - 7y7r) + /% + /%> (5) 

where R^- is the Ricci tensor associated with £(£). In this paper, we raise and lower indices 
on 3-dimensional objects with 7^ and 7^ respectively. 

The ADM coordinate conditions which generalize the isotropic Schwarzschild metric read 

7T U = 0, (6) 

where, and from here onwards, repeated covariant or contravariant indices imply the usage 
of Einstein summation convention and 

7iJ =^% + ^ T , (7) 
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where ip is a conformal scalar and hf? the transverse-traceless (TT) part of the 3-metric 
7^ with respect to the euclidean 3-metric 6ij. By definition, hj^ satisfies hj^ = hj^ = 0. 
The restriction hj^ = gives a simple expression for the 3-dimensional curvature scalar, 
7 1 / 2 R = SipAip, where A stands for the 3-dimensional euclidean Laplacian. The differen- 
tial equation, used for gauge fixing 7^, reads 

3 7yj - ln,i = °- (8) 

Taking into account the gauge condition for 7T y ', namely Eq. (0), the following decomposition 
can be achieved, 

7r« = 7f« + tt^ t , (9) 
where 7f JJ denotes the longitudinal part of 7T W . It may be expressed as 

^=<i + <i-^ m m, (10) 

which implies ir %3 ■ = 7f y ^ = A7r* + |vr^-, suggesting that Eq. (jSJ) can be used to compute 
an elliptic equation for ic l . The TT-part ti 1 \ namely j^- the canonical conjugate 

to /i^ T , which gives the independent field degrees of freedom. The Hamiltonian, which 
generates the time evolution of the independent degrees of freedom of the system (matter 
plus gravitational field), is given by ^| 

E[MV, hff, 4r] = -^Jd 3 x A4,[MV, 4tI (H) 

where MV denotes the (non-specified) matter variables. If the matter system consists of 
black holes, the matter variables MV in our paper enter via boundary conditions at the 
apparent horizons and at spacelike infinity. 

The Eqs. (JHJ) and (J2J) result in the covariant trace of 7r y of the form tt = TT^hJ^ . Taking 
into account the space-asymptotic properties tt 1 ^ ~ 1/r 2 and hj^ ~ 1/r, the gauge condition 
Eq. (6) turns out to mean asymptotic maximal slicing. The gauge conditions Eqs. (6) and 
(7), or (8), are very close to the well-known Dirac gauge conditions. The condition Eq. (8), 
e.g., is identical with the corresponding Dirac gauge condition to linear order in 7^ — 5^. 

The functions ip and 7r l , and hence n l \ are determined using the Hamiltonian and mo- 
mentum constraints, given by Eqs. (j2J and (J3J), by elliptic equations. The elliptic equations 
for the (scalar) lapse a and the (vector) shift functions (3 l result from the Eqs. (J3J) and (JHJ), 
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applying the coordinate conditions Eqs. © and (J7J) respectively. The explicit Poisson-type 
equations for a and (3 l can be derived from 

7 « | m - «' = - a(H - 7 R) - 2«7 vr vr m - 27 7 vr # , m 
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+ a (Y l Y m T lm - 7 il a 2 T 00 ), (12) 

where use has been made of Eq. (J2J . The elliptic equation for (3 l results from 

(/%),; + (%),j - g(%),i = -(«7 _1/2 (2vr ij - 7^-Tr)),,- + i(a 7 - 1/2 (27r, j - 7,'^)),*. (13) 

The functions, hJ 3 T and 7r^ T , in general, follow from the evolution equations, given by Eqs. 
(HJ) and (0). Alternatively, /i^ T and 7r^ T may be obtained from the Hamiltonian, Eq. ([lip, 
via the corresponding Hamilton equations. 

III. THE MINIMAL NO-RADIATION APPROXIMATION 

Motivated by the MHD approximation in electrodynamics, a minimal truncation of the 
Einstein theory, in view of suppressing radiation emission, is achieved by putting only the 
TT-part for the evolution equation for n l \ given by Eq. to zero. The TT-part of Eq. 
(jlj) may be denoted as 

^TT,0 = ^TT ) (14) 

where Alf T = A*- 2 — A^. The expression for A l \ obtainable from the right hand side of Eq. 
after using the Hamiltonian constraint equation, reads 

A ij = - a-f 1/2 (R ij - j lj R) - a 7 - 1/2 (27r im 7r^ - tttt") + 7 1/2 (a |ij - 7^a |m |m ) 

+ (7r^/? m )| m - - 7T™^' |m + ^ a 7 1/2 (7 li 7 jm T im - 7 ^ 2 T 00 ). (15) 

The longitudinal part of A %: > , namely A l [ , may be written in terms of a vector V 1 as 

^£ = V j ,i + ^4 - \^V m m (16) 

which satisfies 

a\" • a«. (17) 

The elliptic equation that determines hJ T , derivable using the condition Aji T = 0, reads 



2 7 V2 

27% = - — —7^7^ 



A In I „ln am\ < „lm on < „nm nl 
A L~V K P )|m + 7T P| m + VT P ]n 



+ 7y(27T n vr m - 7T ) 



2 7/„ _.|m \ 1 .j r . 2-HHK 



-A^jm + 2 7T7Tjj + -l(a w - 7^' | J + — 7(?y + H^T™) . (18) 



The left hand side of the above equation, which introduces the Laplacian for h]j, takes the 
form 



27% = -x 8 
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-Jkj,il - Jil,kj) + ^ll kl l nP (Xn.il^p.kj - ^n.ijFp.kl) , (19) 

where x 8 = ^ ~ l^m^mn anc ^ where the Christoffel symbols Ti j k are given by \ (7^^ + 
lik,j — ljk,i)- For the derivation of Eq. (fTH|) . we have also used the following relation 



77 y — X 8 $ij ~ ^hjj 1 + hJk T hIJ 14 • m terms of a function 0, vanishing at spacelike infinity, 
ip may be given by (1 + |). 

In the Hamiltonian formulation, the present approximation to Einstein's equations is 
identical with the replacement of only the evolution equation 



which is the TT-part of the Eq. (0J), through 



0. (21) 



The dropping of vrrf T0 exactly corresponds to the dropping of the displacement current, 
- in Maxwell-Ampere equation of electrodynamics. This is so, since — E is the momen- 
tum conjugate in the Hamiltonian formulation of electrodynamics. However, this is not quite 

ri 

the way MHD like approximation to general relativity is implemented in p], as they impose 
hJ? Q = in their scheme. In the post-Newtonian approximation, their scheme will not gen- 
erate the correct leading order expression for 7ivj? T . Moreover, the dynamics associated with 
their approximation will coincide with the Einsteinian one through second post-Newtonian 
order only. 

It is obvious that stationary spacetimes are reproduced exactly in our approximation. The 
constraint and the evolution equations, given by Eqs. (j2J), (jHJ), ©, and © still determine 
if), 7i\ a, and (3 l respectively through elliptic equations. Moreover, Eq. (fTS|) defines hj?, 
which describes the independent field degrees of freedom, through an elliptic equation, (see 

n 

Eq. (3.8a) in 4], for more details). In addition, the longitudinal part of ti 13 , defined via V\ 
can be determined through an elliptic equation for V 1 , given by Eq. (17). The remaining 
quantities n^ T should be determined, solely algebraically, using the evolution equation for 



7ij, given by Eq. (|5|). in the following manner. Let us write the Eq. (5) in the form 

7ii,o = % , (22) 

where 

% = « 7 - 1/2 (27r ij - 7lj 7r) + /% + (23) 
Using the definition of 7^, given by Eq. (7), we obtain, 

^ T o = B tJ - -BuS^ (24) 

where B a = 3( 7 mn 7 mn ) -1 (.B" — 7 d 7 4fc ^ T ) with given by — a7~ 1/ ' 2 7 K 7r+2/? J l i . Using these 
equations, 7r^ T follows in terms of hj? , the only time derivative appearing in the minimal 
no-radiation approximation, given by 

0,1/2 

tt^ = - r j - - — c ijlk \f3 llk - /^"VV^VT 1 ] 

a 

+ tlc^ij 1 ^ - 7 ir 7 fcr 7 sm 7 sn (7 M 7 w )->^,o, (25) 
where C^ lk is the inverse matrix of the matrix 

= \ (SuS jh + S lk S 3l - ( 7 ij - 7 l "7 Jri 7 mm (7 P V ,? )- 1 )^ T ) • (26) 

Notice, the matrix C~^ k deviates from the unit matrix in the quadratic order of hJ k T only. 
For maximal slicing condition, = r y v Kij = n-y~ 1 / 2 /2, the resulting C^j k would be a unit 
matrix, when applied to symmetric tensors. 

The internal consistency of the proposed approach comes from the fact that only the 
evolution equations for the true field degrees of freedom are affected. This implies, in the 
evolution equations for the TT variables, 

7T TT0 — ^ ^TT> {^(&) 



c 4 m 



■" lo7rG S-k^t 

if the first equation is solved for hj^ using Eq. (JTHJ), then the second equation can be 
solved for 7ivj? T in terms of hj^, as given by Eq. (|25[). However, putting hj^ = in Eq. 
(I27b|) means destroying the Legendre transformation, which is a fundamental property of 
the physical theory, whereas setting 7r^ T = in Eq. (|27a|) implies a dynamical condition for 



the extremum (minimum) of energy. Imposing hJ 3 T = 0, instead of solving for it, results in 
the well-known Wilson-Mathews approach, if 7ivj? T is determined through Eq. f)22p. We also 
note that a recent effort, the so-called CFC+ approach, which tries to improve conformal-flat 
approximations uses Eqs (|18j) and (J25|) . restricted to the leading post-Newtonian order, to 
get hj? and . This approach is identical with the full Einstein theory through the 

second post-Newtonian order. 

In general, the Hamiltonian for the matter system is the Routh functional depending on 
the matter variables only 

U[MV, t^(MV), K™(MV)\ =H-^ J ^s4^, (28) 

In the stationary case, TZ and H coincide. For our quasi- stationary approximation, we also 
shall take 1Z as Hamiltonian. It agrees with the non-dissipative part of the Einsteinian one 

n 

through third post-Newtonian order as e.g. given in Ref. 

Finally, we point out that hj? is the only partial time derivative appearing in the coupled 
elliptic equations, present in the minimal no-radiation approximation. While employing any 
iterative procedure to solve our coupled elliptic equations, it should be natural to equate 
the (in general) small quantity hj? to zero to obtain the first order solution to hj^ and 
the rest of the variables involved. During the second stage of iteration, one should then 
employ hj?, computed using the first order hj?, to get the second order solution to hj? and 
other quantities. In this manner, the iterative procedure may be extended to higher stages. 
The above procedure was employed, in the post-Newtonian framework, first to compute hj^ 
to the second post-Newtonian order, which was employed to compute hj^, required at the 



third post-Newtonian order calculations 
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This communication will be followed by a detailed article 



where the explicit expres- 



sions, in terms of the basic variables and their spatial derivatives, for the elliptic equations 
for if), a, /3\ and hj^ as well as for ir l and V 1 will be presented. 
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